arXiv:1502.02229v2 [cond-mat.soft] 18 May 2015 


Motility-induced phase separation and coarsening in 

active matter 


Giuseppe Gonnella a , Davide Marenduzzo b , Antonio Suma c , Adriano Tiribocchi d 

a Dipartimento di Fisica, Universita di Bari and 
INFN, Sezione di Bari, via Amendola 173, Bari, 1-70126, Italy 
b SUPA, School of Physics, University of Edinburgh - Edinburgh EH9 3JZ, UK 
C SISSA - Scuola Internazionale Superiore di Studi Avanzati, 

Via Bonomea 265, 34136 Trieste Italy 
d Dipartimento di Fisica e Astronomia, Universita di Padova 
via Marzolo 8, Bari, 1-70126, Italy 


Abstract 

Active systems, or active matter, are self-driven systems which live, or func¬ 
tion, far from equilibrium - a paradigmatic example which we focus on here is 
provided by a suspension of self-motile particles. Active systems are far from 
equilibrium because their microscopic constituents constantly consume energy 
from the environment in order to do work, for instance to propel themselves. 
The nonequilibrium nature of active matter leads to a variety of non-trivial in¬ 
triguing phenomena. An important one which has recently been the subject of 
intense interest among biological and soft matter physicists is that of the so-called 
"motility-induced phase separation", whereby self-propelled particles accumulate 
into clusters in the absence of any explicit attractive interactions between them. 
Here we review the physics of motility-induced phase separation, and discuss 
this phenomenon within the framework of the classic physics of phase separation 
and coarsening. We also discuss theories for bacterial colonies where coarsening 
may be arrested. Most of this work will focus on the case of run-and-tumble 
and active Brownian particles in the absence of solvent-mediated hydrodynamic 
interactions - we will briefly discuss at the end their role, which is not currently 
fully understood in this context. 
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1. Introduction 

Active systems are self-driven systems which live, or function, far from ther¬ 
modynamic equilibrium. A paradigmatic example is provided by a suspension 
of interacting self-motile particles in equilibrium with a thermal bath at tem¬ 
perature T. This system, like active matter in general, is characterized by the 
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continuous conversion of internal energy into work or movement 0, 0, which 
drives it far from equilibrium even in steady state [dj. 

Nature offers many examples of active matter, at very different scales: these 
range, for instance, from the cytoskeleton of eukaryotic cells, to whole bacterial 
colonies and algae suspensions, and even to bird flocks and schools 0, i i & 
0 i, i 111 [Til l . Self-propelled units can also be artificially realized in the lab in 
many different ways, for example by anisotropic surface treatment of colloidal 


particles [12|, ll3|, |14|, US, IlS HH, M 


The inherently nonequilibrium nature of active matter endows it with a num¬ 
ber of non-trivial features that have no analogue in passive, equilibrium materi¬ 
als, such as a suspension of Brownian colloidal particles. For example, bacterial 
swimmers may accumulate near walls due to activity alone fiol. fio| . and their 
large scale diffusive motion may be rectified by using funnels or other asymmet¬ 


ric geometrical obstacles [21, 22|. Within a bacterial suspension, it is also well 


known both that there are violations of the fluctuation-dissipation theorem [23 


and that the diffusion of passive tracers is dramatically affected by the presence 


of active swimmers 



The existence of a parameter to 
for active matter has been anal- 


be interpreted as an effective temperature [3__, 
ysed in |3l], |H, HI, [li], [34]. The interplay between activity and the rodlike shape 
of bacteria and most other active particles leads to further surprising phenom¬ 
ena, such as the large scale coherent turbulent-like motion observed in bacterial 
monolayers (35|, 24, 36, H0 H3, S, Idol ]- or the strikin g l y n on-Newtonian rheology 


found in simulations of model active gels [41], [42, 43], |44], I45| . 


Our main focus in this review is another example of a phenomenon that is 
unique to active systems, the so-called “motility-induced phase separation 1 11 , which 
has attracted a lot of attention recently 06], 03, 03, 49, 50, 51, 03, H3, 54] . This 
term refers to the phase separation into a concentrated and a gas-like dilute phase 
observed in a system of self-propelled particles interacting solely via steric or ex¬ 
cluded volume repulsion, when the overall density is large enough 0. The physics 
of motility-induced phase separation (henceforth MIPS) can be understood as fol¬ 
lows. Imagine that due to a fluctuation the local density of self-propelled particle 
increases in some part of the system. The particles in that region will therefore 
slow down due to the enhanced crowding. Now, as we will understand more 
formally in Section 2, self-motile particles (such as bacteria) accumulate where 
they move more slowly, much as pedestrians in a busy street (and very much 
unlike colloidal particle undergoing Brownian motion). This potentially triggers 
a positive feedback loop, whereby particles accumulate where they are slower, 
further slow down due to the crowding, then accumulate even more etc. The 
theory which we present in Section 2, and which is fully confirmed by Brownian 
dynamics simulations of self-propelled spherical particles (discussed in Section 3), 
predicts that these clusters should coarsen so as to lead to complete phase sep- 


1 We will not consider in this review models, like the Vicsek model (l| , where self-propelled 

particles interact via explicit alignement interaction and phase transition is related to sponta¬ 
neous symmetry breaking - see also the end of Section 3.3. 
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aration. Indeed, the theory suggests that this nonequilibrium phase separation 
should be very similar (in the same universality class) to the phase separation 
in a system of passive colloidal particles which are subject to a mutual attrac¬ 
tive interaction - according to this view the main role of activity is to provide 
an effective attraction between particles, in the absence of any thermodynamic 
attractive force. Indeed the macroscopic description of MIPS is very similar to 
a Cahn-Hilliard equation for a phase separating binary fluid, albeit with some 
subtle important differences. 

Recently, experiments have also been reported where MIPS was observed 


in suspensions of self-motile particles only subject to steric interactions m 

m- 


Clustering of self-propelled synthetic particles was also found in Ref. |56l. |57 
although in both those cases the coarsening appeared to arrest, for reasons which 
are still not fully understood. In particular, it is not clear to date to what extent 
the theory of MIPS (as discussed above and in Section 2) applies to the work of 


Ref. 15(1 1571 158| . or how to modify it in an experimentally relevant way so as to 


lead to microphase separation, or arrested coarsening. 

Another important aspect of MIPS worth discussing is that the theory nor¬ 
mally considers spherical self-propelled particles, whereas in most cases active 
particles are in practice elongated (think for instance of bacteria or chemically 
powered nanorods). It is therefore of interest how MIPS appears in a suspensions 
of rodlike, or elongated, rather than spherical, microswimmers. While several 


simulations of rodlike active particle exists |59l. l60t l61| . it is only recently that 
the effect of particle shape on MIPS has been directly addressed - for instance 
with simulations of dumbbells [52j. We will review these works in Section 3. Sec¬ 
tion 3 also contains a brief description of coarse grained theories for self-propelled 
particles which align locally (for instance due to steric interactions if they are 
rodlike in shape); such theories are more complicated that the Cahn-Hilliard-like 
equation governing the evolution of active spherical particles because they need 
to follow the dynamics of the polarisation held as well. 

Finally, we will consider other ways in which MIPS can be extended in Sec¬ 
tion 4. First, we will discuss an intriguing and generic way to arrest phase 
separation, by coupling the motility-induced feedback described above to a lo¬ 
gistic law for the population of active particles, used as a simple metaphor for 
bacterial reproduction |62], The arrested patterns predicted by the resulting 
theory closely resemble those found experimentally in growing bacterial colonies 
within semi-solid media such as agarose gel, which are normally explained in¬ 
stead by bacterial chemotaxis 63|. Second, all theories discussed so far neglect 
solvent-mediated hydrodynamic interactions between active particles. However, 
in a nonequilibrium active suspension, hydrodynamic can potentially modify 
the steady state attained by the system, and not only its transient dynamics. 
Furthermore, unlike active particles in Brownian simulations, neutrally buoyant 
hydrodynamic swimmers in a solvent need to be force free, so that the minimal 
active contribution is a force dipole (sj. In Section 4.2 we will briefly review some 
of the studies which have appeared in the literature for suspensions of spherical 
squirmers 64, 65, 166!, Isif, where the activity arises through a non-zero imposed 
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slip velocity at their boundary. The consequence of solvent-mediated hydrody¬ 
namics for MIPS is currently unclear, and some conflicting conclusions have been 
reached in the literature, possibly reflecting the dependence on details of near 
held interactions. 


2. The physics of motility: self-trapping of self-propelled particles 


In this Section we provide a simple theoretical framework to understand the 
theory of MIPS, or self-trapping of self-propelled particles BET 


50, 21 . In 


particular, we will describe this phenomenon in terms of an instability appearing 
in the coarse-grained “hydrodynamic” equation manifested by the presence of a 
negative diffusion coefficient. To establish a hydrodynamic equation for the self- 
propelled swimmers, we use here a kinetic approach; similar calculations have 


been also presented in Ref. [10| and are an adaptation of Refs. 68, 1691 


2.1. A coarse-grained theory of motility-induced phase separation 

Let us consider a population of swimmers that runs for a time r with velocity v 
between two successive tumbles (sudden changes of directions [Tl]). Let f(x,u, t ) 
be the distribution of active particles with position x, swimming in the direction 
u at the time t. The kinetic (Boltzmann-like) equation for the evolution of / in 
d dimensions can be written as 


d t f(x,u,t) 


V x ■ [vuf(x,u,t )] 


f(x,u,t) J f{x,u,t)d d l u 

T T J 

solid angle f 1d 


( 2 . 1 ) 


The first term on the right hand side of the equation is a convective term, where 
vu is the velocity; note that in general we are interested in the case in which v 
depends on the spatial position, x. The second term is a linear term accounting 
for the particle that "tumble out" of u , with a rate 1/r ; the same rate appears 
in the third term, that accounts for the particles "tumbling into" the direction 
u. The integral of / appearing in the third term is the zero-th order moment 
of / with respect to the orientation u , which is nothing but the density p of 
swimmers in x. Finally, denotes the d- dimensional solid angle. 

A way to analyze Eq. (12.Hi is to consider the so-called hydrodynamic limit 
in which suitable moments of the distribution functions, which are linked to 
macroscopic observables, change smoothly in space and in time. The idea is to 
perform a gradient expansion (which is very similar to the so-called Chapman- 
Enskog expansion (T0|) of Eq. (12.11) . which assumes that successive derivatives 
(gradients) are smaller and smaller on macroscopic (coarse grained scales). The 
distribution functions is then written as 


/ = / (0) + / (1) + / (2) + ... 


( 2 . 2 ) 
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where f (n> is a contribution which is ~ (D(D n ), where D denotes derivatives 
over space or time. We assume that in the hydrodynamic limit every successive 
derivative applied to / makes it much smaller, that is: 

C> ( D n ) <C £> (D n ~ l ) < ... < O (D 1 ) < O (D°) = I, 

hence 


f {n) ^ 1) ^ ^ /(l) < /(0). (2.3) 

Then, we can set up a hierarchy of equations at various orders in D, by equating 
terms of order D n in the left and right hand side of Eq. (12.11) . as follows: 


O(D 0 ) : 

f( o) _ P 

1 V d 


(2.4) 

0(D l ) : 

d t f {0) + V* (vufW) = 

/ (1) 

r 

(2.5) 

0{D 2 ) : 

+ V, (vufM) = 

/« 

r 

(2.6) 


From Eq. (12.41) we find that no term other than contributes to the coarse¬ 
grained hydrodynamic density p. Then, by summing Eqs. (12.51) and (12.61) . and 
taking the integral over the directions we get 

d t p + J dp (vupfW) d d ~'u + J dp (vupfW) d d ^u = 0, (2.7) 

where the relations J f^(x, u, f)d d_1 w = f f^ 2 \x 1 u 1 t)d d ~ 1 u = 0 have been 
used (Greek indices denote Cartesian components, and the usual convention of 
summation over repeated indices is implied). Then we can use the expressions 
of /(°) from Eq. (12.41) and fd) from Eq. (12.51) in order to obtain 

dtP+ TT d J 7 M /3P) d d ~ lu ~ 


The second and the third integral in the above expression are zero by symmetry 
so that, using 

I UaUpd^U 5 a0 

^ ■ 7’ (2 - 9) 

we are left with 

d t p = • [vV(pv )]. (2.10) 


r 

fh/ 


r 


d t d a ( vu a p ) d d 1 u 


— y (vu a dp ( vupp )) d“ 1 u = 0. (2.8) 
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Eq. (12.101) is at the basis of the physics of self-trapping, and has some re¬ 
markable consequences H The important point here is that the running velocity 
can depend on the position, that is v = v(x). In steady state conditions, the 
distribution (obtained by setting dtp = 0), is p = l/v(x) - hence bacteria, or self- 
propelled particles, accumulate where they go slower. This may seem intuitive, 
but is an important property which has no counterpart in passive suspensions of 
Brownian particles. Indeed, consider a system of non-interacting Brownian par¬ 
ticles, which have a position dependent diffusivity, D(x). Because the system is 
passive, in equilibrium the probability distribution (hence the density) needs to 
be proportional to the Boltzmann weight, and this would give a uniform density 
in absence of any interaction. Mathematically, this is because the correct coarse 
grained equation for the density of colloidal particles p(x) would be 

dtp = V • [D{x)Vp ], (2.11) 


i.e. the diffusion and the velocity enter in different positions with respect to the 
gradients! As a result, in steady state the only possible solution is p equal to a 
constant, at variance with the case of the self-propelled particles. 

It turns out that this seemingly technical difference leads to dramatic dif¬ 
ferences in the physics of motile and Brownian systems. Consider for instance 
the case of a concentrated suspension of active particles, where the swim speed 
decreases for density, for example, simply due to crowding. Let us consider for 
simplicity an exponential decay with density, 


v ~ e xp/2 v 0 , 


( 2 . 12 ) 


where A measures how steep is the decay of the swim speed with p. Eq. (12.101) 
in ID then becomes: 


dtp = d x 




Or [£Wp)9rf>] ■ 


(2.13) 


The effective diffusion coefficient D e ff is now a function of the density. More 
importantly, for p high enough D e ff < 0, so the effective diffusivity becomes 
negative! Negative diffusion means that an infinitesimally small fluctuation leads 
to a divergent increase in the density held. In other words, the system phase 
separates. A regularised version of Eq. (12.131) can be written down by introducing 
a fourth-order derivative, as follows, 


d,p = d x lD elf (p)dzp ] - kd*p 


(2.14) 


The latter term is analogous to the surface tension between particle-rich and 
particle-depleted domains appearing in standard continuum theories for multi¬ 
phase fluids [71). 

Physically, the phase separation described by Eq. (12.141) can be understood as 
due to the positive feedback mechanism described in the introduction. Due to the 


2 We note that the theory we have derived implicitly assumes that there is no macroscopic 
polarisation of the self-propelled particle orientation. 
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fact that the swimmers accumulate where they go slower, the local density will 
increase leading to more decrease in velocity and eventually to phase separation. 
In reality this happens for a sufficiently steep decrease of swim speed with density, 
as predicted by Eqs. (12.13H2.14j) . This mechanism, based on the appearance of an 
instability threshold for the diffusivity, has been named “self-trapping”, and the 
phase separation itself is sometimes called activity-induced, or motility-induced, 
phase separation (3, 461. 

Eq. (12.101) and Eq. (12.141) have been derived for run-and-tumble particles, such 
as bacteria; another important class of self-propelled colloids is that of “active 
Brownian particles” (henceforth ABPs). In ABPs, the direction of swimming 
only changes gradually, by rotational diffusion (as opposed as to “instantaneous” 
tumbles) - indeed the term Brownian comes from the fact that the typical ori¬ 
gin of the gradual change of swimming direction is rotational diffusion. Like 
run-and-tumble particles, ABPs also exhibit motility-induced phase separation 
(see Section 3, and, for instance, the theory and simulations in 29|). It is 
important to note here that the same theory leading to Eq. (12.101) and Eq. (12.141) 
can be modified in a rather straightforward way to apply to ABPs, essentially 
by substituting 1/r with (d — 1 )Dr [ 46l fij ] (Dr is the rotational diffusion coef¬ 
ficient). 


2.2. Relation with Cahn-Hilliard equation, noise and interfacial terms 

The approach adopted so far neglects hydrodynamic noise, and it further 
lacks a systematic derivation of the —kd^p term, or indeed of any interfacial 
terms, which should be of fourth-order in gradients. 


Here we briefly sketch the theory, initially derived in Ref. 46, 69| and after¬ 


wards generalized in Ref. [501.l72fl . in which these aspects have been included. 

First, we note that Eq. (12.101) can also be written formally as a Cahn-Hilliard 
equation, 


d t p = -V ■ \-M(p)Vp] 


(2.15) 


In Eq. (I2.15j) . the chemical potential p can be obtained from the functional 

differentiation of an effective free energy T with respect to p, with F — f fodr, 

and with / 0 = p(lnp — 1) + In [v(p')]dp' (this was first noted in [46]). The term 

2 

M(p) is an effective mobility, equal to for run-and-tumble particles (recall 
from Section 2.1 that r is the run duration, equivalently the inverse tumbling 
rate). An explicit calculation of the chemical potential as p = SF/Sp and a bit 
of algebra show that Eq. (12.15ft is the same as Eq. (12.101) . 

Within the effective Cahn-Hilliard description in Eq. (12.151) . the role of fluctu¬ 
ations can be included by added a conserved noise as follows (see j46j for details), 


3 tP = -V • \-M(p)Vp + v / 2M(p)A 


(2.16) 


On the right hand side of the equation A is a noise vector with components 
(Aj(r, t)Aj(r', t')) = 8 ij 8 (r - r')5(t - t'). 
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Importantly, a full analysis of the phase separation kinetics requires the pres¬ 
ence of interface-like terms to stabilise domain walls between the two phases. 
This is achieved by assuming that v is not strictly local but v = v(p), with 
p( r) = p + 7 2 V 2 p. This is the leading order correction allowed by rotational 
invariance (50[. Here 7 is a smoothing length proportional to the run length 
for run-and-tumble particles, or to the persistence length of ABPs: in formulas 
one expects q(p) = 70 T r v(p), where 70 is of order unity and 77 is the orien¬ 
tational relaxation time (equal to either r or to D^ 1 ). This assumption re¬ 
lies on the fact that a single ABP samples the density p over a length scale 
much greater than the interparticle spacing and is proportional to the persis¬ 
tence length l(p) = T r v(p)/(d — 1). 

By using the nonlocal form of v(p) = v(p) just advocated, it can be shown 
that the chemical potential including interfacial terms may be written in the 
following form |50j 


p — ln p + In v(p) — K,(p)W 2 p + 0(V 4 p), 


(2.17) 


where n(p) = — 


70 Tr 

(d~ 1) 


<p)%- 


It is important to note that the interfacial term 


breaks the mapping to a “standard” Cahn-Hilliard equation, in that there exist no 
free energy functional of which Eq. (I2.17p is a functional derivative. Equivalently, 
the interfacial term breaks detailed balance. 

Finally, we note that in Ref. j5c| a repulsive term of the form f rep = uQ(p — 
Pi)(p—piY is summed to f 0 to model excluded volume interactions, and to prevent 
the existence of a phase of infinite density. The repulsive free energy term is 
phenomenological, and it consists of constant uj multiplied by the Heaviside step 
function 0 ; there is also a density threshold pi at which the repulsion is switched 
on. 


To make calculations in practice with the continuum model in Eq. (12.161) . 
an explicit expression for both v(p) and M(p) is needed. These require some 
information from microscopic theories, or simulations, which we will discuss next. 
Suffice here to say that, as shown in Ref. p5oj, an appropriate choice at low and 
intermediate densities is the following: 


v(p) = v 0 (l - v 0 a s T c p), (2.18) 

for v(p), and 

M ( p ) = = A)P( 1 - v 0 a s T c p) 2 , (2.19) 

for M(p). In the two equations above, a s is a scattering cross-section. This 
quantity is obtained by assuming that a single ABP trajectory consists of straight 
runs with speed v$ punctuated by collisions lasting for a time r c in which the 
particle arrests its motion. (We note that Eq. (l2.18D is also consistent with the 
limit of Eq. (l 2 . 12 l) for A —> 0.) 

I 11 Fig. [T]we show intermediate-time snapshots of the concentration 0 = p/po 
(where po = (v 0 a s r c ) _1 ) obtained by numerically solving Eq. (12. 161) in 2D (left 
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Figure 1: In this figure we report snapshots for the local packing fraction 0 = 
p/po, with p 0 = ( v o cr s T c)~ 1 , obtained by numerically solving Eq. (12.1 6ft in 2D (left 
panel) and in 3D (middle and right panel), both at a value of the overall packing 
fraction 0 O — 0.5. The right panel shows a 2D projection of the corresponding 
3D density field, whose isosurfaces are shown in the middle panel. More details 
about the mapping between 0 and p are described in Ref. [50|. 


panel) and in 3D (middle and right panel) both at overall packing fraction 0o = 
0.5 [50, 111]. Isolated domains of the dense phase in 2D (in yellow), clearly 


distinguished from the diluite phase (in blue), represent regions where particles 
tend to accumulate. Especially in 3D, the domain topology resembles one which 
is expected of phase separation in passive gas-liquid systems, or binary fluids. 

It is quite a remarkable and striking feature that the model in Eq. (12.16(1 . 
although almost fully built within the framework of equilibrium statistical ther¬ 
modynamics H, captures very well the existence and dynamics of motility-induced 
phase separation, which is a far from equilibrium phenomenon. (Indeed, as we 
shall see more in detail in Section 3, Eq. (12.161) also reproduces well the dynam¬ 
ics observed in direct particle simulations of concentrated suspensions of active 
Brownian particles, with an enourmously reduced computational cost ['72].) 


3. Motility-induced phase separation in self-propelled particle systems 

In practice, both in simulations and in experiments, motility-induced phase 
separation has been observed in active brownian particles, rather than in run- 
and-tumble bacteria (recall that in ABPs rotational diffusion replaces tumbling, 
see Section 2). ABPs can be realised in practice, for instance, as synthetic (bi- 
)metallised colloids where some chemical reaction creates a self-phoretic local 
chemical motoi0 jl2l. 13 . ITdl. [lil . 16 . 17 . 181. An example of a motility-induced 
clusters observed experimentally in such systems is shown in Fig. [21 It is also rela¬ 
tively straightforward to set up a Brownian simulation of ABPs, and because this 


3 The weak transgression from equilibrium and detailed balance is the interfacial term, which 
has been extensively discussed in Ref. [HU. This has been shown to have no impact on the 
coarsening dynamics (see Fig. 5). 

4 Self-phoretic motors are particles which move due to the gradients of a chemical which 
they produce themselves. 
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Figure 2: (a) An example of the dynamics of a cluster of active colloids. The 
particle marked by the black arrow (left snapshot) leaves the cluster and is substi¬ 
tuted by another particle approaching (right snapshots). Red arrows denote the 
orientation of the caps, (b) The mechanism used by a particle for self-trapping 
is shown. After a collision, a particle becomes free when its orientation changes 
due to the rotational diffusion. This figure is taken from Ref. [55| where large 
enough Janus particles, allowing to visually resolve the single polarizations, are 
used. 


is the approach which has been most commonly employed in the study of MIPS, 
we focus in this Section on numerical work, reviewing first the model, then the 
phase diagram for MIPS, and later on the dynamics of MIPS in ABPs (Sections 
3.1-3.2 describe simulations for spherical particles, 3.3 and 3.4 describe more 
general models for rodlike particles, or self-propelled particles with an alignment 
interaction). 

In the simplest theoretical framework leading to MIPS, N spherical ABPs 
interact with each other only by excluded volume effects and are subject to a 
propulsive or active force continuously acting along a fixed polar axis for each 
particle H. This force is typically constant in modulus but its direction changes 
with time due to rotational diffusion. A minimal theoretical model with these 


5 In a system with momentum conservation, as discussed later, the total force on a neutrally 
buoyant swimmer should instead be zero. However Brownian dynamics theories, e.g. simu¬ 
lations, neglect fluid-mediated interactions so the only way to propel a particle is to apply a 
force along its direction. 
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characteristics is the following 


r % 

9i 


ViV 0 + H F h' + vj , 




Vi 


(3.1) 

(3.2) 


Here r* are the coordinates of the center of mass of the particle i (i = 1, ...,N) 
and /i is the mobility, while the polarity axis of each particle is defined by 
z>, = (cos 6 t , sin 6j ), where the angle 6 1 * is measured with respect to some fixed 
laboratory axis. Furthermore, w 0 is the fixed self-propulsion velocity and F^- 
is the force due to the particle j acting on the particle i that takes into ac¬ 
count excluded-volume steric effects. TfJ(t) and rji(t) are Gaussian, white noises 
with zero mean and correlations < vlaifyljgit') >= 2 D8 i:j 5 a p8(t - t') and < 
Vi(t)Vj(t') >= 2 Dn8ij5(t — t'), where D = fcgT/i and Dr respectively denote 
the brownian translational diffusion constant (T is the temperature of the bath) 
and the single particle rotational diffusion constant Dr = 3D/a^, with a c the 
diameter of the colloid. 

The occurrence of phase separation in a realization of the above model was 
first shown in jT7!], and later on confirmed and analysed in several works (which 
are reviewed later on). A nice feature of particle-based studies of MIPS, whether 
experimental or numerical, is that these allow a clear identification of the mi¬ 
croscopic mechanism leading to the initial formation of clusters of self-propelled 
particles. In particular, one observes that when two or more particles collide 
head-on they are blocked due to the persistence of their orientations (see the 
experimental pictures in Fig. [2]). A particle can escape from the cluster if, due 
to rotational diffusion, the direction of the polar axis changes so that the ac¬ 
tive force pushes away the particle from the cluster. This microscopic view is 
complementary to the macroscopic theory of MIPS, presented in Section 2. 

Relevant adimensional numbers in terms of which the above model can be 
analyzed are the following. The surface fraction is given by 

0o = iV § (3.3) 

with S c the area occupied by an individual colloid and S the total area of the 
box where the particles move. Then we introduce the Peclet number, Pe, which 
is a dimensionless ratio between the advective transport rate and the diffusive 
transport rate. For particle flow one defines it as Pe = Lv/D, with L a typical 
length, v a typical velocity, and D a typical diffusion constant. We choose L —y 
a c , v —> v 0 , and D —>• D = ksTn). Then, 

Pe = t 

k B T n 

While other definitions of the Peclet number based on rotational diffusion are 
possible, we use here the same definition given in the references discussed in this 
section. 


(3.4) 
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Figure 3: (Color online.) Phase diagram of a system of spherical ABPs (red con¬ 
tinuous curve), and of dumbbells (blue dashed curve). Inside the curve, at high 
Peclet numbers, the system undergoes phase separation into two phases charac¬ 
terized by two different densities. Here, typically, large and stable clusters are 
observed, as shown in the snapshot for spherical ABPs on the right. In the case 
of dumbbells, these are frozen due to steric interactions and point preferentially 
towards the center of the cluster, while active colloids can still freely rotate inside 
the cluster. For small Peclet numbers the system does not show the formation of 
such large and stable clusters; a typical snapshot of the system in this uniform 
phase is also shown, again for spherical ABPs. The location of the transition line 
is based on the results of simulations fixing Vq = 0.1 in both cases and varying 


see 


5ll . |52 ] for further details); the right snapshot is taken at T = 0.01, the 


one on the left at T = 0.05. Spherical ABPs follow an equation of motion with 
the same parametrization of the center of mass of dumbbells. 


3.1. Phase diagram 

To begin our review of particle-based simulations of MIPS using the model 
described above, we discuss determinations of the phase diagram. From the 
considerations made above, and in Section 2, we expect that there should be 
a phase transition (MIPS) at low temperature or large self-propulsion (i.e. at 
large Peclet number) when the density of the system is sufficiently high. The 
phase diagram of a system of self-propelled disks is shown in Fig. [3J Here, the 
physical meaning of the red curve is similar to that of the coexistence curve in 
fluids with liquid-vapour transition. Outside the coexisting line, the system can 
only be in a single phase, while for Peclet numbers higher than a critical value 
and densities inside the coexistence curve, it phase separates into a gas phase 
and a particle-rich phase. An approximate analytical estimate for the phase 
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diagram can be found in Ref. [60] - while the theory in Section 2 predicts phase 
separation, it disregards thermal diffusion so strictly works only for very large 
(actually infinite) Peclet number. 

An interesting quantitative prediction of the phase diagram in Fig. [3] is that 
the critical point occurs for a value of the Peclet around 50, which can quite 
readily be realised with e.g. bacterial swimmers (the challenge though is to 
concentrate them up to 0 O ~ 0-6, such that MIPS is predicted to occur in Fig. ED- 

To determine numerically the coexistence curve one may, for instance, con¬ 
sider the probability distribution function P of the local density q i>. This can 
be obtained dividing the full system in cells with linear size larger (typically 10 
times) than that of each particle and calculating the density in each cell. The 
size of the cells has to be much smaller than the linear size of the full sample and 
big enough to sample (f) correctly. Then, letting the system evolve from an initial 
homogeneous configuration at a given total surface fraction, one can measure the 
distributions of local densities, P(<f>), in steady state. 

An example of the analysis of these distribution functions, for a system of self- 
propelled disks is shown in Fig. 21 (the corresponding parameters are given in the 
caption). In each panel of this figure one observes that at low Peclet number the 
distribution is characterized by a single peak approximatively corresponding to 
the value ot the total density while at high Peclet the distribution is characterized 
by the presence of two peaks which can be taken as the density values for the 
two coexisting phases. This can be shown to be a consistent procedure, e.g. by 
checking that starting from a different total surface fraction the values of the two 
peaks remain the same 52, [T3J . 




Figure 4: Local density distributions for spherical ABPs at Peclet numbers given 
in the keys in the form Pe(T) for different temperatures, u 0 = 0.01,0.1, cr c = 1, 
and /i = 0.05. The global density of the system is <f >o = 0.5. 


3.2. Kinetics of phase separation 

We now turn to the review of the studies of the kinetics of MIPS (these are 
less numerous in the literature). In analogy with the study of domain growth 


13 





and coarsening in gas-liquid systems and binary fluids, the idea here is to mea¬ 
sure the size of the phase separating domains (the particle-rich ones) L(t), for 
instance by taking appropriate moments of the structure factor [71] . The size of 
these domains then increase in time, and the analogy with Cahn-Hilliard models 
suggests that such a dynamics follows a power law, 


L(t) ~ t° 


(3.5) 


where a is the exponent which measures how quickly coarsening proceeds. It is 
useful to recall here that the passive counterpart of the ABP model we are con¬ 
sidering is a gas-liquid system in the absence of solvent-mediated hydrodynamic 
interactions - under these conditions one expects a = 1/3 (7l| . 


In their ABP simulations, Redner et al. 


49j, measured a growth of the cluster 
mean-size compatible with an exponent a = 0.272. Stenhammar et al 
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a variant of the model of (13.1113. 21) . performing more extensive simulations with 
5x 10 5 particles, found a = 0.28 in 2 D systems. The same authors also performed 
3D simulations with 4 x 10' particles obtaining in this case a = 0.34. 

The previously mentioned power-law regimes were found after an initial tran¬ 
sient with faster growth. It was observed in [72) that the cross-over to the final 
regime occurs when L(t) exceeded a persistence length related to the characteris¬ 
tic balistic displacement of a single particle. This corresponds to the persistence 
length mentioned in Sect.2.2 and can be evaluated as l = vq/(Dr(cI— 1)). Since 
this quantity is proportional to the Peclet number (it depends on the ratio Vq/T) 
and phase separation occurs with Pe > Pe c , it follows that very big systems are 
needed in order to have a large enough extension of the late time regime with 
reliable measurements. For example, in the simulations of [j2j|, Pe=300 implies 
l = 50cr so that a linear size Lbox = 350a was considered leading to 4 x 10' 
particles in order to have 0 O = 0.5. 

In general, the results in both 2D and 3D are interpreted as consistent with 
the a = 1/3 exponent, which is expected for very long times (such times are 
reached in 3D, but not in 2D due to the reasons just discussed). A natural 
question is whether the run-and-tumble particles previously discussed behave in 
the same way. In this area, the largest simulation to date is that in Ref. ji74| . 
which uses a lattice model of run-and-tumble particles and demonstrates that 
the “diffusive” exponent a = 1/3 is found in simulations of MIPS in 2D. 


3.3. MIPS with non-spherical particles: active dumbbells 

Until now, we have focussed our review on the case of spherical self-propelled 
particles. In active matter physics, however a spherical shape as assumed in 
the previous section, is the exception rather than the rule. Swimmers, whether 
synthetic or natural, are often rod-like or elongated: this is true both of most 
bacteria, but also for chemically propelled nanorods jil, 3, @, 75] . It is therefore of 
interest to analyse how MIPS is affected when the self-propelled particles have an 
elongated shape, i.e. they are rods, or dumbbells. In particular we will consider 
in this Section the case of active dumbbells. 
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Figure 5: In this figure the length-scale L(t) for the 2D (left snapshot) and the 
3D (right snapshot) systems are shown. Both of these are obtained for packing- 
fraction 0o — 0.5. The black curves (top set of curves, starting from the leftmost 
point in the figures) refer to L(t) as found by particle based simulations (ABPs); 
the red curves (middle set of curves in the figures) refer to continuum model 
simulations according to Eq.(2.16); the blue curves (bottom set of curves in 
the figures) refer to L(t) for the case in which the detailed balance is restored 
(see |54j|). The dashed lines gives the corresponding fitted exponents, which 
are: a 2 D = 0.279 and a 3 n = 0.341 for ABPs, a 2 D = 0.287, a 3 £, = 0.333 for 
the continuum model with detailed balance violation. r r is the orientational 
relaxation time defined in Sect. 2.2. The dotted lines represent the persistence 
length which separate the region of superdiffusive behaviour (a > 1/3) at short 
times from the region of diffusive behaviour (a ~ 1/3). Figure is taken from 
Ref. |72|. with permission. 


In Ref. |29] a model of active dumbbells was introduced to describe the dif¬ 
fusive experimental behaviour of a bacterial bath coupled to colloidal tracers. 
Other simulations of rodlike active particles were presented in [159, 60, [§Ij: all 
these cases demonstrate that the physics of anisotropic active particle is richer 
than that of their isometric counterparts, for instance due to the possible onset 
of local orientational (nematic) ordering. To the best of our knowledge, the first 
analysis of the effect of shape anisotropy on MIPS was performed in |52j. In 
that work, the phase behaviour of active dumbbells of variable overall density 
and Peclet number was studied by means of Brownian dynamics simulations. 
The model of active dumbbells is a simple generalisation of those of ABPs; how¬ 
ever, for dumbbells the rotational diffusion is no longer explicitely described as 
in Eq. (13.21) . rather it effectively results from the net balance of forces on the 
two disks of each dumbbell (a detailed description of the model is given in [34]). 
Rotational and translational diffusion properties of active dumbbells have been 
studied in 00 3- 

The elongated shape of dumbbells leads to some important differences with 
the case of spherical active particles. To begin with, phase separation is overall 
favoured in the case of dumbbells 0 ]. For instance, in Fig. [3] the coexistence 
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curve for a dumbbell system can be compared with that of colloids; while the 
curves cross close to the critical point, for sufficiently large Pe the density range 
whether MIPS can be observed is larger for dumbbells. 

However, by far the most striking feature about the dumbbell aggregates is 
that they exhibit coherent and long-lived vortex-like rotation (see Fig. [6]). This 
coherent rotation is absent in the motility-induced clusters formed with spherical 
particles. On the other hand, we also observe that the clusters observed in the 
system of active rods of Ref. [77l . [78 1 60], where the length-to-width aspect ratio 


is larger compared with that of our dumbbells, show a translational motion but 
little detectable rotation. 

What drives the self-organisation of rotating aggregates in systems of dumb¬ 
bells? And why are they absent for spherical swimmers? A first explanation for 
the formation of such structures comes from considering torque balance on the 
droplets which form during MIPS. Within the Brownian dynamics simulations 
performed in j52|, the forces exerted by each dumbbell add a local torque which is 
balanced by the drag on the cluster. While also spherical self-propelled particles 
would exert a local torque, this can only lead to rotational diffusion and not to 
sustained rotations, because the local swimming direction of each active particle 
can rotate freely within the cluster [[49]. This is not the case for the dumbbell 
aggregates, where steric interactions essentially quench the polarisation of each 
of the dumbbells in the rotating droplet (see Fig. ED- 

The overall coherent rotation of the cluster in Fig. [H] is apparent from a coarse 
grained plot of the dumbbell momentum held pv (see Fig. [6] B), where p and v 
are the local density and velocity held of the dumbbells at a given point in 
space and instant in time. It can be seen that the topology of pv is that of a 
vortex. Inspection of the local ordering within the high density rotating clusters 
of dumbbells also demonstrates that there is both nematic and hexatic order 52 


The local director prohle, measuring the average orientation of the dumbbells, 
P. shows a characteristic spiral pattern (Fig. [6p). A distinguishing aspect of 
the rotating aggregates in dumbbell systems is therefore that the velocity and 
orientational patterns are distinct: the former is a vortex, the latter a spiral. A 
minimal mean held framework to understand this fact was proposed in Ref. (52| : 
a conclusion from this theory is that a requirement in order for polarisation and 
velocity held not to be parallel to each other is that the theory should include a 
non-trivial momentum balance between self-propulsion and steric repulsion. 

Computer simulations suggest that the velocity of the rotating clusters, fl, 
should depend on their radius, R. In particular, the Brownian dynamics simula¬ 
tions of Ref. ||52| measure a scaling of H with 1 /R. This scaling can be understood 
by means of dimensional analysis as follows. First, we note that the torque on 
the cluster must be proportional to R 3 F a , where F a is the magnitude of the ac¬ 
tive force acting on each dumbbell, as the individual dumbbell torque is ~ RF a , 
and given that the polarisation within the aggregate is coherent (interestingly, 
other arrangements, e.g. random, would lead to different scaling). Second, the 
rotational friction scales as R 4 |52j. Balancing the two contributions leads to the 
observed 1/R scaling. It is interesting to highlight here that rotating clusters 
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Figure 6: (A) Snapshot of a phase-separated active dumbbell fluid with three 
clusters (two of which are touching). The clusters, from the smallest to the largest 
one, contain about 4000, 12000, and 25000 dumbbells, respectively. The inset 
shows a detail of the dumbbell configuration; red and green beads indicate the 
tail and head of each dumbbell and are connected by a line. (B) Coarse grained 
velocity held corresponding to the snapshot in (A), obtained by local averages 
on a square mesh. (C) Coarse-grained orientation profile of the dumbbell fluid 
(P), obtained by averaging as in (B), again correspondin g to the configuration 
in (A). (D) Hexatic order parameter (computed as in Ref. j49j j) corresponding to 
the configuration in (A) (defects, or grain boundaries, betweeen ordered regions 
are observable). Figure is taken from Ref. [52:]. 


have been found in experiments with mixtures of bacteria and polymers, where 
bacteria feel an effective attractive interaction due to depletion attractions aris¬ 


ing from the presence of the polymer [79(. In the experiments, again a 1/R 


behaviour of the angular velocity of the clusters is observed, however the details 
of the balance leading to this scaling are different. In Ref. [79], bacteria aggre¬ 
gate in 3D, the torque comes from bacteria at the surface of the cluster, and the 
drag on the cluster is Stokesian and ~ R 3 . Furthermore, the rotating clusters 
found in Fig. [H] appear due to MIPS, whereas those in |79| due to an additional 
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interparticle attraction (the depletion attraction). 

While the work in Ref. [ : 52j covers the effect of shape anisotropy on the MIPS 
phase diagram and cluster morphology, an interesting open question is to what 
extent and how shape may affect the kinetics of phase separation, e.g. the coars¬ 
ening and growth laws of the domains. 

Finally, it can be observed that active dumbbells differ from spherical ABPs 
because their shape introduces an additional hydrodynamic variable, the local 
polarisation P. In this sense, they can be viewed as a particular example of 
a more general class of self-propelled particles, which interact with each other 
through an alignment interaction. In the case of dumbbells, this interaction 
comes from steric effects alone, but in other cases it may have different ori¬ 
gins: for instance bacterial swimmers may interact due to hydrodynamic inter¬ 
actions [3|. In the best-known model for active particles, introduced by Vicsek 
et ah in 1995 111, an ensemble of self-propelled particles interact solely via an 
alignment interaction whose origin is unspecified. The Vicsek model exhibits a 
transition from an isotropic phase to a flocking state with active particles moving 


coherently (80), |8l|, 82|, |83|, 84|. A review on different dynamical features in the 


Vicsek model is beyond the scope of our current work and the interested reader 
can see e.g. Ref. (85|. A generalization of the Vicsek model in which effective 
steric interactions among active particles are incorporated is studied in Ref. [86s |. 


4. Other aspects of motility-induced phase separations 

4.1. MIPS with reproducing particles: arrest of coarsening 

In Sections 2 and 3, we saw that the physics of motility-induced phase sepa¬ 
ration is akin to that of phase separation in a gas-liquid system, or in a binary 
mixture. This is because the clusters, which form when active particles bump 
into one another, coarsen, e.g. by Ostwald ripening - i.e. smaller aggregates 
evaporate to increase the size of larger clusters. In some cases, however, MIPS 
can also be arrested, so that clusters coarsen until a certain equilibrium size is 
reached, and no more. Here we consider one interesting example in which this 
happens: a population of reproducin g r un-and-tumble self-propelled particles 
whose velocity decreases with density |62|. 

A simple way to describe self-replication in bacterial colonies is to do so via a 
logistic growth law. This can be done naturally within our hydrodynamic frame¬ 
work, and leads to the following mean-held equation for the bacterial density, 


d t p = d x [D eff (p)d x p\ - kd*p + ap (1 - (4.1) 

where a is the growth rate (in s" 1 ), and po is an equilibrium density which the 
bacteria tend to in the absence of other interactions. The value of a may vary 
greatly within bacterial colonies, between 0.6 —0.8 x 10 -4 s" 1 for rapidly growing 
bacteria, to 0.1 — 1 hr" 1 , or even less for slowly growing ones. The slowdown 
of swim speed with density can be driven either by quorum sensing between 
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the microorganisms, or by steric effects. The former is probably more likely in 
bacterial suspensions, as typically the equilibrium density that they would grow 
to, po ; which sets a scale for typical values of p, is below 1%, at which it is 
unlikely that steric interactions play a major role. 

Eq. (14.ip leads to an arrest of the phase separation, due to the additional re¬ 
production (and death) time scale which is set up by the logistic term, ap ^1 — -j^J 
The mechanism leading to pattern formation and selection is as follows. First, 
MIPS leads to the formation of clusters, or equivalently to an inverse diffusive 
flux which drives bacteria, or active particles, from dilute regions into particle- 
rich ones. This flux is counteracted by the logistic term, which tends to decrease 
the local density if it is above the target value po (indeed p = po is the only sta¬ 
ble solution for bacteria with a swim speed which is constant, or which decreases 
only mildly with density |62|). The competition between these two terms leads to 
a steady state where the bacteria accumulate only to a certain point, and where 
multiple clusters coexist in steady state. One important finite length scale, which 
is to a good approximation related to the intercluster mean separation in steady 
state, is given by a Jvjoi, where T> is an effective diffusion constant (depending 
on the density imbalance between the particle-rich and the dilute phase). Only 
if a —» 0 (no reproduction), coarsening can proceed indefinitely, corresponding 
to macroscopic phase separation discussed in Sections 2 and 3. 

The resulting pattern depends strongly on the initial condition. For instance, 
starting from a uniform phase with noise leads to the formation of droplets, 
whereas starting from a bacterial droplet, or “inoculum” (i.e. with the density 
concentrated in a small region initially) leads to the formation of concentric rings 
or dot patterns, not unlike those which are usually associated with chemotactic 
bacteria (63] (see Fig. d). 

Interestingly, there are other experiments which suggest that, even with syn¬ 
thetic self-propelled particles, motility-induced phase separation results in sep¬ 


arate clusters rather than in macroscopic phase separation |56|, 157], [58|. The 


mechanism leading to arrest of coarsening under those conditions is still unclear. 
Clearly, it cannot be the one discussed previously, because while the overall 
density of synthetic self-propelled particles is conserved in the experiments of 
Ref. 156. 57, IH|, it is not in Eq. (14.Ill , or in the bacterial pattern experiments it 
relates to. 


f.2. The role of solvent-mediated, hydrodynamic interactions 

Up until now we have described models of “dry” active matter [§], i.e. where 
the solvent where active colloids typically move is neglected. In reality, micro¬ 
scopic swimmers (which move at low Reynolds numbers) stir the fluid they are 
in, hence they interact each other also through the solvent hydrodynamic flow 
held. This important fact, not considered in the previous Sections, can have rel¬ 
evant effects on the phase behavior. We will now briefly discuss this point in this 
Section. Models for active matter under fixed how prohle have been considered 


m 
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Figure 7: Patterns formed in a colony of self-reproducing bacteria with density 
dependent swim speed; reproduction and death are modelled by a logistic growth 
law. The initial condition was an inoculum (small bacterial droplet at the centre 
of the simulation box). Each row provides a series of snapshots observed for 
increasing time; the color refers to bacterial density (see colorbar on the right). 
More details on parameter values etc. are found in Ref. [62], Figure taken 
from (62], with permission. 


There are different swimming mechanisms used by microorganisms, artificial 
active colloids and active droplets. Here we consider the most studied case of 
the so-called "squirmers" that propel themselves by a prescribed axisymetric 
surface velocity held (9cj, 91]. They can be used as a simplified generic model 
for diffusiophoretic particles flfil . 92l |. or for biological microswimmers where the 
action of cilia or flagella can be described at mesoscopic scales by the squirmer 


surface velocity held |90l . l9l| . A review on the role of hydrodynamics in the 


self-propulsion mechanisms is given in |93|. Here, we will focus on the effects of 
the hydrodynamic interactions on the collective behavior of squirmers and, in 
particular, on the possibility of having a (motility-induced) phase separation as 
in the cases discussed in the previous sections. 

We hrst show in Fig. [S] the velocity held of different kinds of single squirm¬ 
ers. In order to describe different real microswimmers, the characteristic of the 
model can be varied to describe either pullers (equivalently contractile swim¬ 
mers), with a vortex behind the body (Fig. [8](c) ,(f)), or pushers (equivalently 
extensile swimmers) with the vortex in front of the body (Fig. [H](a) ,(d)). Fig. [S] 
also shows the velocity held for the case of neutral swimmers. The velocity held 
associated with a moving puller or pusher is dipolar (the velocity of the solvent 
decays with distance from the swimmer, r, as r~ 2 ), whereas that associated with 
a neutral swimmer is quadrupolar (the velocity decays as r~ 3 ). 

The simplest way to analyze the hydrodynamic interaction between squirmers 
is to consider the how held around a pair of squirmers kept at hxed position 94 


In the case of pullers, the fluid is pressed into the region between the two squirm- 
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ers increasing the pressure there, and leading to repulsion. On the other hand, 
in the case of pushers, the squirmers are expected to attract each other. An 
analytic expression for the force between two squirmers is available decaying 



ented in parallel, allowed to swim freely, are observed to move towards each other 
but, due to diffusion, they do not always collide [94j. Sometimes, after collision, 
rotational diffusion reorients the squirmers and the trajectories diverge, with an 
increasing average spreading with decreasing Peclet number. 


Pusher Puller 



Figure 8: Flow streamlines generated from an isolated squirmer in the swimming 
frame (top row) and in the lab framce (bottom row). Left and right panels ((a)- 
(d) and (c)-(f)) show a pusher and a puller, respectively, and the corresponding 
dipolar velocity held. The middle panel ((b)-(e)) shows the_quadrupolar velocity 
held of a neutral swimmer. The hgure is taken form Ref. |64l |. 


There are few studies of the collective behavior of squirmers. Simulations 
based on Stokesian dynamics [96] , 95] reported in (64| showed that pullers are 
systematically more ordered than pusher spherical suspensions. In this work, 
order was measured in terms of the global polarization value coming from the 
average of the orientation vectors of all particles. In some way, this result contra¬ 
dicts the expectation based on the attractive (repulsive) character of the inter¬ 
action between pushers (pullers). This behavior was attributed to a prevalence 
of head-to-head collisions favouring a faster decorrelation of pusher swimmers. 
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The result was also shown to be independent on the initial, aligned or isotropic 
conhguration of the ensemble of particles. It also contrasts a prediction from 
continuum theories saying that an isotropic state is unstable only for nonspher- 
ical pushers The number of particles considered in the simulations of (64 ] 
was, however, too small in order to properly discuss the possibility of a phase 
transition. 

A mixed simulation method, with a Lattice Boltzmann scheme for the fluid 
coupled to particle dynamics for the squirmers was considered in [65]. Here the 
number of particles was much larger, few thousands, and the system was 3D. 
It was shown that pullers have a tendency to form transient aggregates. This 
was related to the presence of polar order, which occurs for pullers but not in 
the case of pushers. However, a macroscopic motility-induced phase separation 
was not observed. The origin of the polar order was attributed to a mechanism 
similar to that discussed in ji3] • The aggregation phenomenon was shown to be 
almost absent in the case of pushers. 

Using simulations methods similar to those of |64j], it was shown in [&6[ that 
hydrodynamics should suppress phase separation for isometric squirmers in 2D; 
this was explained by mapping the model to a system of active particle in a 
range of parameters where motility-induced phase separation does not occur. In 
this work, most of the simulations were done with neutral swimmers (although 
simulations suggested that the results were similar for pushers or pullers). 

The same squirmer system was considered in another recent work by Zottl 
and Stark [67], that presents the results of simulations based on the multiparticle 
collisional dynamics method coupled to particle dynamics for the squirmers. The 
swimmers move between two walls, which allow both 2D translations and full 3D 
rotations. Unlike Ref. [66;] , this study found a pronounced difference in behaviour 
in the two cases. The main result of this work is that at sufficient high area 
fraction </> 0 > 0.5 the system can separate in a gaslike and a crystalline phase, 
depending on the nature of the swimmers. Pusher are observed to form a single 
cluster at 0o ~ 0.65 while pullers are observed to form several exagonal structures 
and a single cluster only at very high density. On the other hand, aggregation 
in clusters is seen to be greatly favoured for neutral swimmers and diminishes 
when the pusher character increases. These features are explained in terms of a 
strong enhancement of the rotational diffusion coefficient, compared with that of 
a single swimmer, due to hydrodynamic swimmmer-swimmer and swimmer-wall 
interactions. 

Finally, we mention the work in Ref. |98|, which used a fluid-particle dynamics 
approach, and studied the effect of hydrodynamics on the collective dynamics of 
dumbbell swimmers with prescribed dipolar force pairs in 3D. This work showed 
that within semidilute suspensions of swimmers (in the presence of thermal fluc¬ 
tuations) hydrodynamic interactions appear to enhance the dynamic clustering 
(which can be seen as a precursor of fully-developed MIPS) at relatively small 
volume fractions. Within this model, this result arises due to hydrodynamic 
trapping of one swimmer by another, induced by the active forces. The final pre¬ 
diction of this model regarding MIPS is still pending. Preliminary simulations 


22 



within this framework suggest that motility-induced phase separation should oc¬ 
cur for sufficiently large propulsion forces (at a volume fractions of 0.2 — 0.25); 
however the structure formed was not compact and exhibited large fluctuations, 
possibly because of the relatively small system size considered. A final assessment 
of the fate of MIPS in this model therefore requires even larger scale simulations. 

As is clear from the above, gaining an ultimate understanding of the role 
of hydrodynamics interactions in motility-induced phase separation is a highly 
non-trivial, and still largely open question. This is due both to the fact that near 
field effects, which are difficult to capture accurately, seem to matter, and to the 
need to run very large scale simulations (as in the case of ABP |50l. lo4l| 1 to come 
to a firm conclusion on phase separation in the bulk. In this sense, large scale 
simulations in the future may hold the key to finding the answer to this difficult 
question. 

5. Conclusions 

In conclusion, we have reviewed here the physics and dynamics of motility- 
induced phase separation, a nonequilibrium phenomenon which occurs in a sus¬ 
pension of self-propelled active particles. Motility-induced phase separation was 
first predicted theoretically in [46| (a nice review on the topic, complementary 
to our work here, is the recent Ref. |S|). The physics underlying this phe¬ 
nomenon can be understood as follows. Imagine that due to a fluctuation the 
density of the suspension increases locally: because self-propelled particles (un¬ 
like Brownian colloidal particles) accumulate where they are slow (see Section 
2 and Refs, therein), then this fluctuation can trigger a positive feedback loop 
whereby particles accumulate further, slow down even more etc, eventually re¬ 
sulting in macroscopic phase separation, at least for suitable values of the density 
and propulsion speed (or Peclet number). Microscopically, the formation of dense 
regions can be traced down to the self-trapping of active particles, which move 
much faster than their diffusion, hence bump into each other to form clusters 
from which they cannot easily escape. 

The physics and dynamics of motility-induced phase separation shares sev¬ 
eral aspects with that of gas-liquid systems, or binary mixtures. Importantly, 
the baseline model predicts that the particle-rich domains should coarsen, with a 
growth law which is compatible with that of a diffusive Cahn-Hilliard equation, 
both in 2D and 3D. This conclusion is gained, for instance, from an analysis 
of very large scale Brownian dynamics simulations of spherical active Brown¬ 
ian particles, and is also confirmed in 2D lattice simulations of run-and-tumble 
particles. We have also seen that the phenomenology of motility-induced phase 
separation has some important differences when considering elongated particles, 
such as dumbbells. Here, particle-based simulations show that the aggregates 
which form due to the self-trapping mechanism rotate coherently, and we have 
seen how the dependence of angular velocity on cluster size can be understood 
with simple scaling arguments. 

Importantly, coarsening can be arrested; an intriguing example is provided by 
a population of reproducing run-and-tumble particles with a density-dependent 
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swim speed, which would lead, in the absence of reproduction, to macroscopic 
phase separation. The presence of an additional time scale, linked to reproduc¬ 
tion and death of bacteria, which can be, for instance, captured by a simple 
logistic law, allows phase separation to be arrested. The resulting patterns are 
very similar to those found in bacterial colonies growing in semisolid media, such 
as agarose gel. Such patterns were normally attributed to chemotaxis; the results 
in Section 4 suggest that arrested motility-induced phase separation provides an¬ 
other viable mechanism for these clusters. 

There are several open questions which are needed to fully understand motility- 
induced phase separation, and we have mentioned some of these in our work. 
First, it is unclear to what extent the diffusive scaling exhibited by spherical 
self-propelled particles, and the associated continuum theory, applies to suspen¬ 
sions of anisometric active particles, such as dumbbells or rods. Second, a few of 
the experiments studying the phase behaviour of self-propelled particles (57L ll00l | 
seem to be consistent with arrested phase separation. The possible mechanism 
leading to arrest of coarsening in this context is unclear at the moment. Third, 
the role of solvent-mediated hydrodynamic interactions are currently not fully 
understood; yet the simulations performed to date suggest that they are likely 
to affect the physics of motility-induced phase separation significantly. We hope 
that our review will stimulate further work on this interesting nonequilibrium 
phase transition in active matter. 
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